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Electron Precipitation Pulsations 
% 

F. V, Coroniti and C. F. Kennel 
University of California, Los Angeles 

ABSTRACT 

When high frequency wave turbulence is present, low frequency 
micropulsations can strongly modulate the high frequency wave 
amplitudes, leading to finite amplitude pulsations in the loss 
rate of energetic electrons from the magnetosphere, micropulsation 
amplitude is small. An extremely idealized model suggests that 
the precipitation modulation depends exponentially on the 
micropulsation amplitude, when the micropulsation period is less 
than the electron precipitation lifetime. 
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I . . Introduct i on 

It has boon popular to relate the loss of charged particles from the 
magnetosphere to the atmosphere to pitch angle diffusion driven by microscopic 
plasma -turbulence,. (Brice, 1963; Andronov and Trakhtcngchrts, 1964; Kennel 
and Petschek, 1966; Cornwall, 1966; Roberts, 1968, 1969; Kennel, 1969.) 

Those theories have generally described temporally quasi-steady precipitation. 
However, observed electron precipitation, which is well documented, is 
rarely even quasi-steady, exhibiting a variety of temporal features on 
millisecond to perhaps thousand second time scales (Anderson, 1968). Less 
is known about the temporal structure of proton precipitation. In this paper 
we attempt a plausible explanation for some of the electron precipitation 
modulations which is within the framework of pitch angle diffusion theories. 

In so doing, we propose a new, but simple-minded, nonlinear interaction 

• * 

between wave modes which is stronger than those conventionally discussed 
in the plasma turbulence literature. 

• The shortest millisecond to second time scale fluctuations in the 
• electron precipitation rate may well be due to ionospheric effects near the 
detector (Lampton, 1967; Perkins, 1968). The similarity between the spatial 
morphology and structure of x-ray microbursts and whistler mode chorus 
emissions, pointed out by Russell, et al. (1969), suggests that 1 second structure< 
precipitation is related to the structure in the basic turbulent spectrum 
responsible for precipitation. The mechanism driving precipitation pulsations 
of periods longer than a few seconds need not be localized near the earth, 
and may not depend strongly on the structure of the turbulent spectrum, since 
electrons bounce many times through the equatorial plane in a pulsation 
period, during which time the chorus structure has changed. Hopefully, long 



period pulsations can be explained by modulations of averaged pitch angle- 

diffusion processes, which are localized to the equatorial plane. 

* ‘ 

Quasi-periodic electron precipitation pulsations have a wide range of 
time scales. Various pulsations with periods from 5 to 300 seconds have 
been observed by Anger, et al . (1963), Evans (1963), Barcus, et al. (1965), 
Barcus , et al. (1968), Paths , et al. (1968a), and reviewed by Brown (1966) 
and Anderson (1968). Riometer and balloon x-ray techniques typically 
measure precipitation pulsations of electrons with greater than 30 keV 
energy. However, auroral light also is quasi-periodically modulated, 
suggested that precipitation fluxes of lower energy (1-10 keV) electrons 
pulsate also (Belon, et al. , 1969). The > 30 keV measurements 

indicate that precipitation pulsations only occur on an already enhanced 
precipitation background (Parks et al., 1968a), a fact of crucial importance. 
In the pulsations, the peak- to- valley precipitation intensity ratios are 
typically 1.5-3, so that pulsations are not small perturbations upon the 
background precipitation rate, another important fact. Finally, the peaks 
tend to have a harder energy spectrum than the valleys. 

Using ground-based magnetometers, McPherron, et al. (1968) found that 
at least those precipitation pulsations with periods in the 5-40 second 
range are accompanied by magnetic micropulsations with approximately the same 
■period as the precipitation pulsations. It is tempting to suggest that many 
low frequency precipitation pulsations are caused by micropulsations. 
(However, Maehlua and O'Brien (196S) have proposed a different mechanism for 
50-200 second period precipitation pulsations .) Without good micropulsation 
information, it is impossible to formulate a precise theoiy of precipitation 
pulsations. Therefore, we can at most outline arguments which indicate that 
small amplitude micropul sat ions could produce large amplitude precipitation 
pulsations. 



-4- 


There are several intuitive difficulties vdth the hypothesised r.dcropulsatio 

source for precipitation pulsations. First, 5-300 second micrepulsations , 

\ t 
acting alone, conserve the electrons' first and second adiabatic invariants. 

Secondly, it is unlikely that micropulsation magnetic., amplitudes be comparable 

with the main magnetic field in space. Furthermore, Parks et al. (1900b) have 

observed low frequency modulations of the energetic electron distributions in 

space, but of small amplitude; any associated micropulsation field would then 

have to be small. Thus , micropulsations appear neither fast enough nor strong 

enough to case precipitation pulsations with 2:1 peak-to-valley intensity 

ratios. However, it is crucial that precipitation pulsations occur only during 


active precipitation events when the general precipitation levels are enhanced. 

In the context of whistler turbulence theory, this means that before micro- 
pulsations can significantly affect precipitation, enhanced energetic electron 
injection must increase the trapped Van Allen electron fluxes above the critical 
flux for whistler instability. Then, the micrbpulsations can affect the 
precipitation rate indirectly, by modifying the 1 whistler instability. 

Consider the quasi-steady enhanced electron precipitation state which 
apparently exists before the onset of electron pulsation activity. The electrons 
are in pitch angle diffusion equilibrium with (for example) whistler turbulence, 
with injection of new electrons balance by pitch angle diffusion to the loss 
cone and subsequent precipitation loss. This equilibrium, a delicate balance 
between sources and sinks of particles and waves, keeps the growth rate of 
whistlers near but slightly above marginal stability. However, virtually any 
small external perturbation of the electrons resonant with whistlers, caused, 
for instance, by micropulsations, changes the whistler growth rate. With the 
whistlers already near marginal stability, a linear change in the growth rate 
leads to an exponential change in the whistler amplitude and therefore in pitch 



angle diffusion and precipitation rate. Conversely, when the whistlers are 
below marginal stability, there is no enhanced precipitation background, and 
micropulsations would not markedly affect precipitation. Thus, we wish to 
describe a nonlinear interaction, through the resonant particle distribution, 
between low and high frequency waves which is exponentially large in the low 
frequency micropuls ation amplitude. This, exponential nonlinearity should 
dominate the algebraic nonlinearitics almost exclusively discussed heretofore 
in the plasma literature. By means of this interaction, the electron's first 
adiabatic invariant is violated, though the micropulsation alone would conserve 
it, and large amplitude precipitation pulsations are possible even when the 
micropulsation amplitude is small. 

In Section 2, we review a simplified model of the electron pitch angle 
diffusion equilibrium; here our understanding, though limited, is reasonably 
secure conceptually. In Section 3, we estimate the change in whistler growth 
rate due to an idealized micropulsation perturbation of the electron distribution, 
in which the micropulsation increases the magnetic field strength seen by 
resonant electrons. Other perturbations would give differences in detail, 
but hopefully not in qualitative behavior. (One complication, which we can 
handle to a limited extent, is that the modulated whistler turbulence itself, 
through enhanced pitch angle diffusion, tends to counteract the micropulsation- 
induced distortion of the whistler-resonant electron distribution, thereby 
reducing the net change in the whistler growth rate.) These crude estimates 
of the growth rate modulation permit similar crude estimates of the change 
in whistler amplitude, from which wc determine the electron precipitation 
pulsation rate, in Section 4, Despite the extreme assumptions made, wo arrive 
at' several physical predictions. In particular, the envelope of precipitation 
modulation depends strongly on the ratio of wave period to average pitch angle 
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diffusion time and upon whether or not the peaks and valleys correspond to the 
strong or w oak pitch angle diffusion rates. 

While some electrons may be precipitated by whistlers some of the time, 
not. all electrons can be precipitated by whistlers all the time. While we 
have phrased the present theory in terms of whistler turbulence theory, our 
discussion is sufficiently loose that other pitch angle diffusion theories 
could be straightforwardly modified to include the interaction with low 
frequency waves . 
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2 • Stead y Pre cipita ti on State 

2.1) Intr oducti on_ 

Here v;e review a grossly simplified picture of a steady state pitch angle 
diffusion process , originally discussed by Kennel and Pctschek (1966), herein- 
after called KP. The same simplifications will be used in describing 
precipitation pulsations. We pay especial attention, in 2,3, to the distinction 
between strong and weak pitch angle diffusion; since the precipitation rate is 
relatively insensitive to the diffusion coefficient in strong diffusion, 
precipitation pulsations will be weak in this limit. 

2 • 2) Diffu sion Equ ation s 

Since whistler waves ought to encounter the greatest number of resonant 
electrons near the equatorial plane, and since they may be damped elsewhere 
(Kennel and Thorne, 1967), the turbulent waves 'are likely to be most intense 
near the equator. We assume that all wave particle interactions occur in a 
quasi -uni form region at the equator. Since the electrons must bounce many 
times through the turbulent region for the diffusion picture to hold, our 
conclusions to follow are limited. to time scales longer than a bounce period. 
Since whistlers propagating parallel to the local magnetic field grow faster 
than oblique whistlers, we assume that parallel waves characterize the spectrum 
at the equator. Here, only electrons which satisfy the cyclotron resonance 
condition, v^ = w-JT/K, ’ interact resonantly with the turbulent distribution, 
to is the wave frequency, K the wave number, f? the electron cyclotron 
frequency, and v a the component of electron velocity parallel to the 
magnetic field. A discussion of the effects of oblique whistlers is contained 
.in Kennel and Pctschek (1969). 
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V/hen the parallel electron velocity v J( is sufficiently large, cyclotron 
resonance interactions produce approximately pure pitch angle diffusion. A 
steady state electron distribution balances pitch angle diffusion to the loss 
cone and subsequent loss to the atmosphere with an injection source S . If 
v is the particle speed and a - cos~* v w /v , the pitch angle, the electron 
distribution function f(a,v) obeys 


jL 3_ 

si no.' 0a 


( D(a) slii a S(a,v) 


( 2 . 1 ) 


outside the loss cone. D(a) is the pitch angle diffusion coefficient, to be 

estimated shortly. Let the equatorial loss cone be at a . Inside the loss 

■ o 

♦ 

cone, ct < a Q , v;e neglect the source, but include an atmospheric sink: 


1 V f 

Ct *00. 


( D(a) a J~ ) - r- = 0 


( 2 . 2 ) 


B 


where the time scale for loss is the quarter bounce time Tg . Since averaged 
over many bounce periods, f is even in v„ , similar equations apply in the 
two hemispheres 0 <_ a _< ir/2, t\/ 2 < a < t> . 

The cyclotron resonance condition specifies only the parallel energy of a 
resonant electron. Therefore, a given wave interacts with electrons of different 
total energy. If the resonant region is slightly inhomogeneous spatially, a 
wave can interact with electrons of different parallel energies. Finally, a 
given electron can interact with different components of the wave spectrum at 
different places. It is reasonable, therefore, that as far as the particles 
are concerned, the effective diffusion coefficient behaves as though the waves 
were quasi-uni formly distributed in K . Henceforth, we will neglect structure 
in the whistler spectrum. For a relatively smooth wave spectrum. 


e^B' ^ -1 

DCci) - ~ 2 ~ 2 ~ x (AKv cosa) = 

me 


cosa 


( 2 . 5 ) 
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where e is the electronic charge in esu, m the electron mass, and c the 
speed of light. AK is the effective widths of the wave spectrum, and B' 
is the wave amplitude integrated over the spectrum. 

Equations (2.1) and (2.2) have been solved by KP for the special case 
that the. injection source is concentrated at flat pitches, S(ct,v) - So (a-u/2) 

When c< > a , 

o * 


f 



-- {log(sina./sina ) * h(a ) } 
o 


( 2 . 4 ) 


where h(ct) is the solution within the loss cone 


h(a) « 



h _ 

Y TO 


where I & , are Bessel functions of imaginary argument. 


(2.5) 


•v The diffusion coefficient involves the wave energy. In principle 
therefore we should account for growth and decay following the individual 
ray paths of all components of the turbulent spectrum. However, the wave 
energy losses probably all occur after the waves have converted out of the 
region of growth. Thus a schematic equation describing the balance of growth 
and loss is 


ft B’ 2 = 2(y-v) B' 2 (2.6) 

where y is a "typical” growth rate (defined below) , and v is a "typical" 
loss rate, which is the order of the group velocity divided by the length 
of the growth region. ' Again, we have averaged over all structure in the 


spectral distribution. 

Y » the "typical" whistler growth rate, is roughly 


Y - ftqA 


(2.7) 
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where rj end A are defined in Kp. Physically, n is roughly the fraction 

of the total electron distribution .near resonance with whistlers. A is a 

Measure of the resonant electron pitch angle anisotropy, l'or a two temperature 

anisotropic Maxwellian, A - T x - T ( /T tf . 

y 

In steady state, 3B 1 '/3t must be zero. To Maintain a steady state, 

waves must e-fold a few times in crossing the equatorial plane to make up 

for losses elsewhere; if y < v , there is no equatorial whistler growth, 

D - 0 , and the source increases the trapped flux (q) unopposed by 

precipitation until y = v . If y » v , waves c- fold rapidly, rapid 

precipitation decreases n , and therefore y , until y approaches v . 

Primarily the number of trapped electrons adjusts until y ~ v , since the 

anisotropy is fixed by the boundary conditions of the diffusion problem. 

2 

Since (2.3) indicates that D(a) depends dominantly upon C , 

(2.6) also describes the time dependence of the diffusion coefficient, 
whereupon , ' - 

~ r 2(y-v) D • (2.8) 


2.3) Limits of Weak and Str ong Diffusi on 

The electron lifetime T. is the total number of electrons outside the 

L 

loss cone divided by the diffusion flux into the loss cone. Using (2.4) 
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Tj nay also be computed by taking the directional flux within the loss cone, 
integrating over the loss coho, dividing by T R to give the loss rate, and 
normalizing to the total trapped flux. From (2.4) and (2.5), we compute the 
ratio of average directional flux within the loss cone J to the averaged 

directional trapped flux, J, f , by first computing the omnidirectional 

2 

fluxes and dividing by the solid angle ratio a. /2 . Thus, 


* - 
[ 3 , <S(v)> 

(>v ) v dv “if- 2 - a 

, o ■ J v o 1 0 . 


J T(>v o ) 


• da’ h(a') 


if/ 2 


[ y^dv — [ da' si net'' 

h (a ) + log 

fsinct 1 ] 

*" ** • * — j 


J v D o Ja 

' 0 L 

isrna. 


V 0 C o 


\ ° J 



a 


or 
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2D T t 


. V 


r * I h (“ 0 5 ♦ l°8C2/e a H = f 1 


a 


( 2 . 11 ) 


where T f is defined in (2.9) and T x , - 21^/a , the minimum lifetime, 

to be discussed shortly. 

In the weak diffusion limit, defined by a o /D o Tj, « 1 , an electron, 

once having diffused to the loss cone, is lost on the next bounce to the 

2 ' ’ 

atmosphere. When /D Q Tg « 1 , h (a Q ) « 1 and is inversely 

proportional to the diffusion coefficient, - 3 /D q . Since increases in 
the whistler amplitude decrease the electron lifetime, the whistler amplitude 
can adjust so that precipitation balances injection. Since J^/J^ « 1 , 

the fluxes in and near the loss cone 'are anisotropic. 

2 

When /D^Tj, - 1 , a particle can random walk across the loss cope 
•before it reaches the atmosphere. Not all the particles in the equatorial 
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plane loss cone will then actually be precipitated. When o '/D^T^ » 1 , 
nearly all loss cone particles remain trapped. In this limit, the steady 
diffusion problem approaches one without sinks, for which the only steady 
solution is an isotropic pitch angle distribution. Since the electrons then, 
are equally probably distributed at all pitch angles, the lifetime is the 
quarter-bounce time divided by the probability that the particle be in the 

loss cone (which is just the ratio of the solid angle outside the loss cone 

2 

to inside), or - 2T^/ ow . This is clearly the smallest possible 
lifetime; at L - 6 , for - 40 keV electrons, it is roughly 200 seconds. 

Once the strong diffusion limit is approached, large increases in the 
diffusion coefficient diminish the precipitation lifetime only slightly. 

Since •Ip/J-p 1 > the fluxes near the loss cone approach isotropy as the 
diffusion coefficient increases. 


2.4) Summary 

In the steady precipitation state preceding a precipitation pulsation 

I . • 

event, the sources and sinks of particles and waves are delicately balanced. 
Pitch angle diffusion to the loss cone removes particles at the rate they are 
injected, at the same time adjusting the electron pitch angle anisotropy and 
intensity so that the average wave growth rate balances the average wave 
loss rate. 

In describing this steady state, we arrived at several idealizations to 
be helpful for the time-dependent discussions to follow. Namely, we shall 
neglect all structure in the wave distribution and pitch angle diffusion 
coefficient, assume that the "typical" growth rate (2.7) adequately describes 
the averaged whistler spectrum, and neglect all details of the propagation 
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and loss of whistler energy, lumping all uncertainty into the 
convective loss rate v . 


effective 
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3. Mo dulat ion of Turbulent Whistler Distribution 

A wave exerts a periodic, force which distorts the particle's 
distributions. In particular, a jaicropulsatlon will modulate the 


velocity 

distribution 


of electrons resonant with whistlers. Over part of the ruicropulsation phase, 
the whistler growth rate should increase. Ordinarily, this small increase is 
not significant; however, when whistler turbulence is already present, small 
increases in the already positive growth rate acting on an already enhanced 
wave distribution can produce significant increases in wave energy density. 
Whistler ray paths will also be modified by the micropulsation, altering 
thereby the effective loss rate of whistler energy. However, it is difficult 


to visualize cases in which growth and loss do not become unbalanced; thus 
turbulent levels and electron pitch angle scattering rates should be modulated. 


Without observations of .micropulsation polarizations, the low frequency 
wave forces are ill-defined. As a qualitative .illustration, we consider 
model micropulsation perturbations which modify the magnetic field strength 
felt by the resonant electrons. For example, the micropulsation could vary 
the magnetic field strength, or transport resonant electrons and waves into 
an increasing magnetic field, or both. In this model, high frequency whistlers 
find themselves in a time varying magnetic field B(t) , 


B(t) - B (1 + b sin w t) 
o o 


(3.1) 


where « is. the micrc-pulsation frequency, and b « 1 . 

; We now attempt to relate temporal changes in the typical whistler growth 
Tate (2.7), 


d toy _ 3 ftn B 3 & n A _3 £n_n 

dt 3 1 3c 3t 


s 


(3.2) 
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to B(t) . The growth rate should be modulated simply because the electron 
cyclotron frequency, the basic scale .frequency for whistler processes, changes 
with time. Secondly, the pitch angle anisotropy might be modified. Thirdly, 
because Jp»| changes, the total electron density and the resonant electron 
density should also change. However, q , their ratio, should not change 
as much from this effect. On the other hand, if the loss rate of high energy 
electrons is modulated due to changes in the turbulence level, q will 
decrease. Thus, enhanced turbulence acts like collisions to relax the 
perturbation in the distribution function. However, it cannot entirely 
counteract the basic perturbation, Many other effects cannot be treated 
Within our approximations. For example, changing |b| alters the energies 
of electrons resonant with whistlers; perhaps destabilizing some waves and 
damping others preferentially. 

In the absence of whistler turbulence, micropulsations would conserve at 
least the first and second electron adiabatic invariants . First invariant 
conservations implies that Tj__ is proportional to B . Assuming that T„ 
does not vaiy strongly, we may estimate 

d in A _ Tj - Hn B _ 1 +_A d An_B 

dt~ " T - T 9t "A dt • ^3.3) 

Since the equilibrium A - ■g- (KP, 1966), A is roughly proportional to An B . 
The modulation of the anisotropy when whistler turbulence is present will 
in fact be smaller than in its absence, since pitch angle diffusion relaxes 
anisotropy changes. 

In steady state, precipitation balances injection, and q is constant. 

How ever, increases in the turbulent whistler level (duo to micropulsations) 
unaccompanied by increases in the injection level, might diminish q through 
enhanced precipitation. However, if the fractional decrease of q in one 
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vave period is snail, this effect will. -'hot completely counteract, increases in 

* 

the growth rate. Let T, be the electron lifetime averaged over a micro - 

L 

* 

pulsation period; when w T. » 1 , n will remain roughly constant over 

O L 

a wave period. The observation that the. peak to valley precipitation rate 

* 

ratio is two, implies T^/2 , so that comparison of the micropulsation 

period with the 'equilibrium electron lifetime establishes tho validity of this 

- . ■ * 

small relaxation limit. Mien w T. >> 1 , v;e may drop 8 £n q/9t ; 

substituting (.3,3) into (3,2) implies E } ~yf " * Integrating 

2 

this equation, using (3.1), and dropping terms of 0(b) , 


Y(t) 


1 + sin w t 
A o 


(3.4) 


where we have assumed the equilibrium growth rate y o equals the equilibrium 

loss rate v . 

o 

it 

When w T, « 1 a full diffusion treatment is needed to find y(t) . 

0 1 # • * 

* 

When « 1 the micropulsation slowly modulates the equilibrium 

electron distribution. For a rough estimate here, we use the equilibrium 
formula for q (derived in KP) to find - . Substituting 

this into (3,2), and integrating. 


Y(t) _ 


o 


B(t) 

B 


1/A D 


D(t) 


(3.5) 


Here, time increases in the diffusion coefficient clearly reduce the modulation 

of y(t) . 

In addition to the above relaxation of the growth rate due to enhanced 
precipitation, there might be other changes in y due to "ripples" in the 
pitch angle distribution of scale Act = /to /D arising from nonmonotonic 


pitch angle solutions to the full time dependent diffusion equation. 
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These higher modes should have a faster time scale than , end their 
precise details should vary spatially on scale length of the electron gyro- 
radius. However , since the observations of precipitation modulation integrate 
over many electron gyro- radii, the "ripples" may be smoothed. 


Ti me D ependence of Diffusion Coeffici ent 

We now solve for D(t) using (2.8) with a time dependent growth rate. 
Given the crude approximations leading to (3.4) and (3.5) it is wise to 
attach significance only to the qualitative behavior obtained in several 
limiting cases. • 

First, we consider the small relaxation limit, where w >> D . 

o 

Substituting (3.4) into (2.8), and integrating, assuming P = D^ at t - 0 


D|tl - D o exp 


2v b 

KfA 


(1 - cos o) t) 


(3.6) 


In the small relaxation limit, the modulation of D is exponentially large 
in b . This case gives the largest precipitation modulation. Since the 
increment of whistler energy created in the micropulsation compression phase 
equals that reabsorbed in the rarefaction phase, D(t) is always greater 
than . Relaxation will gradually decrease the growth rate and damp the 
diffusion coefficient modulations on time scales of 1/D . Presumably, as 
more and more particles are lost, more and more waves will sink below marginal 
stability, and the whistler spectrum will shrink. Should the micropulsation 
also act as a particle injection source, relaxation could conceivably be 
counteracted to some extent. Finally, in strong pitch angle diffusion, the 
small relaxation limit ought to apply for all frequencies, since the enhanced 
loss of electrons will be small. 
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Now we turn to the quasi-static limit where D ~ to . Since a particle 

can diffuse to the loss cone in a wave period, (2.8) must bo solved including 

% 

the changes in I) due to particle losses. Substituting (3.5) into (2.8), 

2 

dropping terms of order b ' , we find 


D(t) - D o + 


2v D 
oo 

1 * C2Y 0 /V : 


2v 

° 

— r - sin w, „t - 
w 

o 


1 - COS W Q t 


to 


b 

A 


(3.7) 


If W 0^ V 0 <<c 1 » then 

D(t) = D o (1 + ~ sin w o t) . ' (3.8) 

Here the modulation of the diffusion coefficient is only linearly propor- 
tional to the micropulsation amplitude, and the whistler amplitude responds 
adiabatically to the micropulsation modulation. 

In summary, the whistler growth rate should be modulated by a low 
frequency micropulsation because the resonant electron distribution is 
modified. Over part of the micropulsation phase, the growth rate should 
increase. When the wave period is much less than the electron lifetime, 
enhanced precipitation does not remove many particles in a wave period. 

Here the growth rate is linearly proportional to the micropulsation amplitude. 
The whistler amplitude is exponentially proportional to the micropulsation 
amplitude. When the wave period is much longer than the lifetime, the 
growth rate modulation is much weaker, since the perturbation is more 
nearly counteracted by enhanced diffusion, and the whistler amplitude 
depends only linearly upon that of the micropulsation. 
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4 . Precipitation Modul ation 

The small loss cone appropriate' to many geophysical situations permits 
a simple expression for precipitation modulation, found without solving the 

n , 

2 

full time dependent diffusion equations. If a Q * « D /w , the pitch 

% 

angle profiles in and near the loss cone approach the steady state profile, 
since particles can random walk across the loss cone in one wave period. 

Thus the ratio of precipitated to trapped fluxes, observed near the earth 
at the mirror points, will be given roughly by (2.12), with D a prescribed 
function of time. This steady state solution is valid only near the loss 
cone. In addition, we can only describe the ratio of trapped to precipitated 
fluxes since the absolute , flux levels can be affected by variations in the 
effective source strength of particles to small pitches by time dependent 
ripples of scale angle V5) /w in the distribution function at large pitches. 
Assuming that the time dependence of D dominates all others v;e find 


yt) 

J T ct) 


~ Tiain 
" yt) 


( 4 . 1 ) 


where T^Ct) is the effective electron lifetime assuming that diffusion 
continues steadily at the instantaneous rate D (t) . 

We may now distinguish four separate precipitation modulation envelopes, 
depending upon the background diffusion rate and the micropulsation amplitude. 
Either the background diffusion rate is- weak or strong; similarly, the 
maximum diffusion rate can be weak or strong ,’ implying the four combinations 
listed in Table 1. 
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TABLH 1 


Case Number 


Minimum Diffusion Rate 


Maximum Diffusion Rate 


1 .) 

Weak 

2 .). 

Weak 


Weak 


Strong 


3.) 


Near Strong 


Near - Str 


ong 


4.) 


Near Strong 


Strong 


The weak-weak case can be calculated easily, since the maximum diffusion 
rate is still weak , so that (t) 1/D (t) . Substituting the zero 

relaxation approximation expression for D (t) , (3.6), into (4.1) we find 


yo y n 

yt) ' t l (o)’ 


2v b 
o 

|w A 
1 o 


(1 - COS 0) t) 

o 


(4.2) 


where 1^(0) is the lifetime at t = 0 . Similarly, the near strong-weak 
and near strong-strong cases can also be computed, since the minimum 
diffusion rate is strong. 


yo , % 

ye “ 1 ' d o t b *** 


2v o b 

l" 0 l A 


(1 - COS W ■ t) 


{4,3) 


. . 4v 


where a /D T n « 1 .. When ~ — t—T b « 1 • there are oscillations about 

o o B A I to | * 


: 2 

the background state J /J. = 1 - ct /D T_ . When 

P t O D 


4v 


A[tt) f 

r\ 1 


; -r b » 1 , 


J- /J-, = 1 for significant portions of the micropulsation period. Here the 
Jr * 

peaks of the oscillation are "clipped." The weak- strong case is mixed, with 
the valleys behaving like (4.2), and the peaks like (4.3), with clipped 
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prof ilcs. A sketch of these four coses is shown in figure 1. 

In the quasi-static limit of eq, (3,8), the precipitation profile should 

A ' 

vary more or less linearly, corresponding to slow, relatively small amplitude 
swells of -the background precipitation rate . 

Assuming that any time dependent ripples in the pitch angle distribution 
spatially average out and, in the zero relaxation limit, that the trapped 
flux is approximately constant over one micropulsation period, (4,2) describes 
•the magnitude of the precipitation flux in the weak-weak case, h’e can then 
renormalize (4.2) to the background precipitation rate, (0) , and taking 
the logarithm of both sides we have 


log 



2v b 
o 


(1 - cos w t) 

k! a 


(4.4) 


Log Jp(t) sinusoidally varies with a positive definite time dependence, 

1 - cos to^t . This behavior qualitatively agrees with observations by 
Parks et al. (1968). 

We now estimate b by taking a typical' peak to valley ratio 
Jp ( it) / dp (0) ^ 2 , a background whistler growth rate V Q ^ 1 rad/s cc , 

A 'v 1/6 , and a micropulsation wave period of 10 seconds, to find b ^ .02 . 
Thus the physical statement that when the whistler noise and electron 
distribution are in diffusion equilibrium, a small distortion of the electron 
.distribution can produce large changes in the whistler amplitude and 
precipitation rate appears self-consistent. A direct comparison with 
observations is difficult. For example, an absolutely reliable ground based 
• magnetometer measurement of micropulsation amplitude would not yield the 
amplitude in space, because for instance, attenuation in the ionosphere and 
at the io7tosphc.rc-naut.ral atmosphere interface is difficult .to estimate. 
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Moreover, the mode structure along the lines of force is not known, so that 

the amplitude at the top of the ionosphere cannot be directly related to 
\ * 

that in the equatorial plane where whistler modulation occurs. Similarly a 
satellite measurement of micropulsations in space would not necessarily 
relate to precipitation pulsations unless correlated modulations of the- VLP 
spectrum were also found. 



We summarize our physical argumehts, comparing them with observations 
where possible. Our arguments are far from rigorous, and suggest only 
that large precipitation pulsations could be caused by small micropul sat ions , 

(1) At disturbed times , enhanced injection increases the trapped electron 
fluxes above the critical flux for precipitation instabilities. We 

• have considered whistler instability as an example. 

(2) A diffusive equilibrium is established, in which whistlers scatter 

electrons into the atmospheric loss cones at the rate they arc 

• « 

. injected, and the whistler growth rate from equatorial resonant electron 
interactions balances wave losses from convective propagation. 

(3) More than one instability is likely to occur at disturbed times, for 
instance, in addition to the high frequency whistler turbulence 
driving particle precipitation, low frequency micropulsations are also 
enhanced at disturbed times (Coroniti et al., 1968). 

(4) The low frequency oscillations perturb the diffusion equilibrium by 
distorting the resonant electron distribution and modifying the 
whistler ray paths. When the frequencies differ greatly, as for 
micropulsations and whistler turbulence, this interaction is physically 
visual iz able.. The micropulsation ..adi.abatica.lly modulates the velocity 
distribution of those electrons resonant with whistlers. This velocity 
space distortion changes the whistler growth rate, which depends upon 
the resonant electron distribution. On the other hand, the increased 
whistler turbulence which follows from the increases in the growth 
rate acts like collisions to relax the micropulsation modulation. 
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However, this relaxation can never completely counteract the micro- 
pulsation perturbation, so the whistler growth rate should always be 
modulated at least somewhat by micropulsations. Only when the whistlers 
are already unstable can small changes in their growth rate lead to 
large changes in their amplitude. Thus precipitation pulsations should 
occur only superposed on an already enhanced precipitation background, 
in agreement with observations by Parks ct al, (1968a). 

(5) If the micropulsation period is much shorter than tlve precipitation 
lifetime, only a small fraction of the electrons can be precipitated 
in a wave period, and the unrelaxed growth rate modulation is large in 
this limit. The modulation of the whistler amplitude is exponentially 
large in the micropulsation, amplitude. Since the minimum lifetime is 
a few hundred seconds on auroral lines of force, large whistler 

modulations ought to be found only in the 3 - nl 300 second period range. 

* 

(6) If the micropulsation period is comparable with electron lifetime, the 
modulations of the whistler intensity should be only linearly proportional 
to the micropulsation amplitude. 

(7) Pulsations of the whistler amplitudes in space produces electron 
precipitation pulsations. McPherron et al. (1968) have correlated 
precipitation pulsations with Pc2 and Pc3 micropulsations. 

(8) Unlike the whistler amplitude modulations, precipitation pulsations 
should have a variety of temporal envelopes depending upon the strength 
of the pitch angle diffusion rate, and upon the micropulsation frequency. 
The largest pulsations should occur when the minimum pitch angle 
diffusion rate is weak, and when the micropulsation period' is much 
shorter than the electron lifetime. Even here, however, there are two 
distinct cases. Either the maximum diffusion rate is weak, in which 
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case the precipitation is exponentially Modulated Kith a smooth profile, 
or the maximum diffusion rate is ^strong, in which case the precipitation 
profile is clipped. During the raicropulsation compression phase the 
precipitation rate exponentiates to its strong diffusion Maximum, 
whereupon it stays constant, until the .raicropulsation decompression 


reduces the diffusion coefficient sufficiently that the precipitation 
rate falls exponentially below strong diffusion to its background weak 
diffusion rate. This cycle repeats. When pitch angle diffusion is 


always strong, increases in the diffusion coefficient do not increase 
the precipitation rate. Here, relaxation of the Micropulsation 

i 

perturbation should be small for all raicropulsation frequencies, and 
v/histler amplitude modulations could be large even at low frequencies. 
However, these modulations would not appear in the precipitation. 

Finally, when the micropulsation period and weak diffusion lifetime 

* . 

are comparable, precipitation pulsations should be weak. 

(9) Several features of the pulsation envelopes described by Parks et al. 
(1968a) may be consistent with the ideas presented here, in particular, 
the sinusoidal time dependence of the logarithm of the observed x-ray 
fluxes. There is a suggestion of a clipped precipitation pulsation in 
their figure 4. 

(10) Since low energy particles have larger fluxes than high energy particles, 
they are more likely to be strongly unstable and in strong diffusion. 

In this case, precipitation modulation would be most pronounced in 
higher energy channels. At the very highest energy channels the 
particles are not sufficiently numerous to be unstable in the first 
place. The high energy particles precipitation might still be weakly 
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rnodulated due to, parasitic effects discussed by Kennel end Pets did; 
(196SJ). These energy considerations are summarised in figure 2. 

(11) Experimental correlations between enhanced whistler (or other high 
frequency) noise and micropulsations are needed to test these ideas. 
Helliwell (1965) lias described quasi -periodic VLF emissions with 
periods similar to those of micropulsations and precipitation 
pulsations, and has already suggested that hydroinagnetie modulation 
of the basic VLF source might be responsible. 
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. ; Figure Captions 

Fig, 1, Possible Precipitation Pulsation Envelopes.. 

The 'logarithm- of the ratio of peak to valley precipitation 
fluxes, log.. Jp/Jy > has a time envelope which depends upon the 
strength of the pitch angle diffusion rate and upon the strength 
of the micropulsation perturbation. The four cases, summarized 
in Table 1, are sketched here. 


Fig. 2. Schematic Energy Dependence of Pulsation Amplitudes. 

If -low energy fluxes are on strong diffusion, they will 

have weak pulsations. Intermediate energy electrons (here 

2 

chosen arbitrarily to lie between 10 and 10 keV) could be 
unstable, but on weak diffusion; these Will have strong 
pulsations. Finally, high energy electrons could have weak 
pulsations due to parasitic effects (Kennel and Petschek, 1969). 
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